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Abstract. We present theoretical and numerical evidence for a random matrix 
theoretical approach to a conjecture about vanishings of quadratic twists of certain 
L- functions. 



In this paper we present some evidence that methods from random matrix the- 
ory can give insight into the frequency of vanishing for quadratic twists of modular 
L-functions. The central question is the following: given a holomorphic newform 
/ with integral coefficients and associated L- function L/(s), for how many funda- 
mental discriminants d with \d\ < x, does Lf(s,Xd)i the L-function twisted by the 
real, primitive, Dirichlct character associated with the discriminant d, vanish at the 
center of the critical strip to order at least 2? 

This question is of particular interest in the case that the L-function is associated 
with an elliptic curve, in light of the conjecture of Birch and Swinnerton-Dyer. This 
case corresponds to weight k = 2. We will focus on this case for most of the paper, 
though we do make some remarks about higher weights (see (26) and below). 

Suppose that E/Q is an elliptic curve with associated L-function 



for Sis > 1. Then, as a consequence of the Taniyama-Shimura conjecture, recently 
solved by Wiles, Taylor, QW], [TW]), and Breuil, Conrad, and Diamond, Le is 
entire and satisfies a functional equation 



where N is the conductor of E and we — ±1 is called the sign of the functional 
equation. Note that we have normalized the coefficients a* so that the functional 
equation of Le relates the values at s and 1 — s. The numbers a* satisfy |a* | < d(n), 
where d(n) is the number of positive divisors of n, and a n — ^/na^ is an integer 
(related to the numbers of points on E mod p for primes p which divide n). Let 

1 Supported in part by a Focused Research Group grant from the NSF. 




oo 



(1) 



71=1 




Typeset by 4^5-TgX 



2 



J.B. CONREY J.P. KEATING M.O. RUBINSTEIN N.C. SNAITH 



d represent a fundamental discriminant and let Xd be the associated quadratic 
character (i.e. Xd(n) = (^), the Kronecker symbol). We assume, for simplicity, 
that (d,N) = 1. Then the twisted L- function is 



(3) L E (s,Xd) = Y^ 

71=1 

This L-function has functional equation 



a* n Xd(n) 



(4) ( j^Tf ) T ( s ) L e( s , Xd) = $e{s, Xd) = w EXd (-N)$(l - s, X d); 

it is actually the L-function of another elliptic curve, namely the quadratic twist of 
E by d. 

The conjecture of Birch and Swinnerton-Dyer predicts that the order of vanishing 
of the L-function of an elliptic curve at the central critical point s = 1/2 is the same 
as the rank of the Mordell-Weil group of the elliptic curve. Thus it would be of 
interest to find an asymptotic formula for 

(5) V E (x):= Yl 1 

|d|<z 
WEXd(-N) = l 
L E (l/2, Xd )=V 

since this would potentially be counting how often the twists of a given elliptic curve 
have rank at least 2. (Note that we have restricted the sum to twists for which 
the sign of the functional equation is +1; these L-functions vanish to order at least 
2 because of the symmetry implied by the functional equation.) Goldfcld [G] has 
predicted that Ve(x) = o(x). More specifically, he predicts that asymptotically, 1/2 
of all twists will have rank and 1/2 of all twists will have rank 1; consequently 
ranks 2 and higher should be infrequent. We will give a more precise conjecture 
about the frequency of twists with ranks at least 2. 

Sarnak has predicted that Ve (x) should be about x 3 ! 4 . His reasoning has to do 
with the formulas of Waldspurger [Wa] , Shimura [Sh] , and Kohnen and Zagier [KZ] 
which relate the value of L£;(l/2, Xd) to the Fourier coefficient of a half-integral 
weight modular form. Roughly, 

(6) L B (l/2, Xd )-« £ c £ (|d|) 2 /Vrf 

where ke depends only on E and where the integers cs(|(i|) are the Fourier co- 
efficients of a half-integral weight form. The Ramanujan conjecture for these co- 
efficients predicts that c_E(|rf|) <C |d| 1//4+e for every e > 0. If cs(d) takes on each 
integer value up to \d\ 1 / i about the same number of times for \d\ < x, then it should 
take the value about x 3 ' 4 times. 

Using random matrix theory, we would like to give a conjecture of the form 

(7) V E {x)~b E x 3 / 4 {\ogx) eE 

for certain constants and e_E. The basic idea is to regard the family 



(8) 



T E + = {L E (s,Xd) ■ WEXd(-N) = +1} 
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as an orthogonal family, in the sense of the families introduced by Katz and Sarnak 
([KS1], [KS2]). More specifically, this family conjecturally has symmetry type + . 
Thus, for example, we believe that the statistics of the low lying zeros of the L- 
functions in this family will match the statistics of eigenvalues near 1 of the matrices 
in SO(2N). 

The point of departure for our conjectures is the work of Keating and Snaith 
[KcSnl] and [KeSn2] (see also [BH] and [CF]) which indicates that the moments 



M E (T,k) = ± ]T L E (l/2, Xd ) h 



\d\<T 

(with T* = \d\<T 1) apparently behave like the moments of the charac- 

teristic polynomials of matrices in SO(2N) where N is of size logT. Precisely, they 
conjecture that 



(9) 

where 

for integer k and 



M E (T, k) ~ 3fc (0 + )a fc ( J B)(logT) fc ( fc - 1 



)/2 



fc-i 



9k (o+) = 2 fe ( fc + 1 )/ 2 n 



2l\ 



(10) 



a k {E) 



fc(fc-l)/2 



(i-^ + A fc + fi + 2E + iV 

V v p) \ v v ) 



p+1 p+1 



This conjecture arises from arithmetical considerations together with the fact from 
random matrix theory that the moments of the characteristic polynomials of matri- 
ces in SO (27V), evaluated at the point 1, averaged over the group can be explicitly 
evaluated. Thus, 



(11) 



M (N,s)= f 
Js 



SO(2N) 
N 

= 2 2Ni 



n 



|det(?7-/)| s dU 

T(N + j - l)T(s + j - 1/2) 



i = ir(j-i/2)r( s + j + jv-i) 

where dU is the Haar measure for SO(27V). The connection with g k is that 
(12) M (N,k)^ g k {0 + )N k ^' 2 

as TV — > oo. Note that the formula for g k (0 + ) can be extended to all real k by 



(13) 



y/G{l + 2k)Y{l + k) 
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where G is Barnes' Double Gamma function. 

Continuing to follow Keating and Snaith ([KeSn2], equations (74) - (81)), we 
observe that knowledge of all the complex moments of characteristic polynomials, 
evaluated at 1, in SO(2N) gives complete information about the density function for 
the distribution of values of the characteristic polynomials at this point. Specifically, 
the latter is the Mellin transform of the former: 

(14) p (jv>) = _L f M (N,s)x- s ds 

ZTTIX J ( c ) 

where (c) denotes the vertical line path from c — ioo to c + ioo; this formula is valid 
for all real x. Note that Po(N, x) dx gives the probability that det(t/ — I) = x for 
an element U of SO(2N). 

For small positive x, the pole of Mo (A, s) at s = — 1/2 determines the dominat- 
ing behavior of Po(N 7 x). In fact, we see that 

:= x-V 2 h(N) 

as x — > + . As A — > oo, 

(16) h(N) ~ 2- 7 / 8 G(1/2)tt- 1 / 4 A 3 / 8 . 

An interpretation of the above is that the probability that an element of SO (2 A) 
has a characteristic polynomial whose value at 1 is A or smaller is 



(17) - / x~ l / 2 h{N) dx = 2X 1 ' 2 h{N). 

Jo 



We apply this reasoning to the values of Le(1/2, Xd)- In particular, by (6) the 
fact that the c,E(|d|) are integers implies that these values are discretized. If, for 
example, it is known that 

(18) L E (l/2, X d)<KE/VWl 
then it follows that L£(l/2,Xd) = 0. Similarly, if 

(19) K E /V\d\ <L E {l/2, Xd ) <4:K E /VWl 



then it must be the case that L E (l/2, Xd) = k e /\/ \d\. 

We assume now that the distribution of values of L E (l/2,Xd) will behave like 
the values of the determinants of random orthogonal matrices with some suitable 
restrictions and use this assumption to conjecture results about the frequency of 
vanishing of L E (l/2,Xd)- The restrictions we have in mind are of an arithmetical 
nature. First of all, we want to include the arithmetical factor ak{E). We expect 
that 



(20) 



M E (T, S )~a s (E)M (N, S ) 
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with N ~ logT. Thus, using 

(21) P E (T,x) = -±- [ M E (T,s)x- s ds 

2mx J (c) 

together with (20), we reiterate the conjecture of [KcSn2] (equation (81)) that 
(22) 

ftp-,,, ~ n ro !^;;™ 3/2) 

= a.^x-^hiN) 

should approximate for small x the probability density function for values of Lb (1/2, \d)- 
Of course, this formula cannot be too accurate as we have already remarked that 
the values of Lb (1/2, Xd) are discretized. The precise nature of this discretization is 
somewhat involved; it involves the constant Ke for which we have explicit formulas, 
but it also involves the coefficients cg:(|d|) whose arithmetic nature is difficult to 
describe. Simplistically, we would like to use the integral of (22), as in (17), and 
(18) to predict that 

(23) #{\d\ <T : L E (l/2, Xd ) = 0,L E (s,Xd) e Fe+} 
~ \^*-i/2T* /T 1/4 h{N) 

with N ~ logT. However, the c(\d\) are divisible by some predetermined powers 
of 2 which change this discretization. For example, in the case of the congruent 
number curve L 32 : y 2 = x 3 — x, the number c(\d\) is divisible by r(d) for squarefree 
d where r(d) is the number of divisors of d. Thus, if 

(24) Ml/2, Xd ) < K E32 T(d) 2 /V\d~l 

then L Es2 (l/2,Xd) = 0. If we introduce the factor r(d) into (23) it will raise the 
expected frequency of vanishing by a factor of about log T giving a total order of 
magnitude T 3//4 (logT) n / 8 for the frequency of vanishing. We expect this to be the 
correct order of magnitude but are not able to say yet what constant we expect. 

If we restricted to prime twists (\d\ — p) then the extra powers of 2 are not so 
significant. This leads to 

Conjecture 1. Let E be an elliptic curve defined over Q. Then there is a constant 
ce > such that 

(25) l~c B T 3 / 4 (logT)- 5 / 8 . 

p<T 
L E (l/2,Xp)=0 

Wc will return to a discussion of the value of ce in another paper. 

We remark that our arguments apply equally well to newforms / of weight 4 with 
integral Fourier coefficients and even functional equation. Here we have a different 
discretization and expect that 

(26) J2 l^c/TV^logT)" 5 / 8 . 

p<T 
L / (l/2, Xp )=0 
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Note in particular that the exponent on T is now 1/4. For a newform of weight 6 or 
higher, we expect that there will be at most a finite number of twists that vanish. 
We have some numerical evidence to support these conjectures. 

For the remainder of this paper we would like to discuss a numerical experiment 
which allows us to skirt the delicate issue of the arithmetic nature of the c(|d|). For 
a prime p we consider the ratios 

(27) R P (T)=( £ l)/( £ i). 

\d\<T \d\<T 
Xd{p) = l Xd(p) = -1 

L E (l/2,Xd)=0 L E (l/2, X d)=0 

L E (s,Xd)^ E + L E (s,Xd)^3 r E + 

By considering this ratio, the powers of T, of logT, and the constants intrinsic to 
the curve E should all cancel out. 
More generally, let 

(28) Q P (fc)=lim L E (l/2, Xd ) k / £ L E (l/2, Xd ) k 

1 — >oo — — 

|d|<T \d\<T 

Xd(p) = l Xd(p) = -1 

LE(s,Xd)eJ r E+ L E (s,Xd)eJ r E+ 

assuming that this limit exists. What is a reasonable conjecture for Q p {k)l 

Using standard techniques from analytic number theory (see [I]), we can evaluate 
Q p (l). Based on the analysis involved in such an evaluation, we expect that 

(29) Q,W = ( ( P+ ' + a i 

(P + 1- a p ) k 

where a p is the p-th Fourier coefficient of the modular form associated with E. The 
heuristics are as follows. Consider either sum that appears in (28) 

(30) £ MVW= E (± a -^ 

\d\<T \d\<T ^n=l 

Xd(p) = ±l Xd(p)=±l 
L E {s,Xd)^^ E + L E (s,Xd)eF E+ 

= ^2 ^ Kxd{n) 

\d\<T n=l H 
Xd(p)=±l 

L E (s,xd)eJ r E+ 



where 



n—nin2---nk 



the sum being over all ways of writing n as a product of k factors. Summing over 
d, the main contribution to (30) comes from those n's that arc of the form p r m 2 , 
i.e. a power of p times a perfect square (since, unless (d, n) > 1, these always have 
Xd{n) = Xd(p) r while, for other n's, we get cancellation as we sum over d). So, the 
main contribution to (30) is roughly 

Xd(p) r 



p r m 2 

\d\<T p r m 2 
Xd(p) = ±l 
L E (s,Xd)£f E+ 
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But, 

b U v — b u b v when (u, v) — 1, so the inner sum above equals 

(i 

Now, 



m z ^ — ' p' 

(m,p) = l r=0 



fc 



Further, one can write an explicit formula for a p r in terms of a p and p. Assume 
that E has good reduction mod p. Writing the corresponding Euler factor (1 — 
a pP s + P 1_2s ) = (1 — a P s )(^ — fip~ s ), with a + P = a p and a[3 = p, we have, 
using partial fractions, a p r = (a r+1 — [3 r+1 )/(a — (3). Substituting this into (31), 
and summing the geometric series we get 



a P-Xd{p) \ 



k oi. 

pZk 



r=0 



P r ) (P + 1 - Xd{p)a p ) 



k ■ 




Hence, the only difference in the numerator and denominator of (28), asymptoti- 
cally, is a factor of 

(p + 1 + a P ) k 
(p+l-a p )k- 

Thus, by the random matrix theory considerations above (as in (15)), taking 
k = —1/2 leads to 

Conjecture 2. With R P (T) defined as above, and E having good reduction modp, 



R p = lim R P {T) -- 

T — >oo 

Note that the number N p of points on the elliptic curve E over the finite field 
F p of p elements can be computed as 

N p = p + 1 - a p 

so that the above ratio is the square-root of the ratio of the number of points on 
E/F p to the number of points on E x /F p for any character \ with x(p) = — 1. 

We conclude with some numerical evidence to support these conjectures. We 
consider three elliptic curves which we call En, Eig, and E 32 - These are associated 
to the unique newforms of weight 2 and levels 11, 19, and 32 respectively. The 
values of L(s,Xd) were evaluated by computing the c(|d|)'s of the corresponding 
weight 3/2 forms. For the level 11 and 19 curves, we only computed these for 
d < and d odd. For the level 32 curve, we computed these for all odd d. The 
relevant form for the level 32 case is described in [Ko]. The forms for the level 11 
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Table 1. A table in support of Conjecture 2, comparing R p v.s. R P (T) 
for the three elliptic curves E n , E w , E 32 (T equal to 333605031, 263273979, 
930584451 respectively). More of this data, for p < 2000, is depicted in 
the figures below. The entries for p = 11 and p = 19 are explained by 
the fact that we are restricting ourselves to twists with even functional 
equation, weXcl{—N) — 1. Hence for En and Ei 9 , we are only looking 
at twists with Xd(H) = Xd(19) = —1. 

and 19 cases can be computed according to [Gr] and were given to us by Fernando 
Rodriguez- Villegas. 
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"E_1 1 R_p/R_p(T) (or p<2000" 



"E_19 R_p/R_p(T) forp<2 



"E_32 R_p'R_p(T) for p<2000" 



Figure 1. Pictures depicting R P /R P {T), 
for p < 2000, T as in Table 1. 
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1e+08 1.5e+08 2e+08 2.5e+08 3e+08 

"E_32, all prime twists" + 



1e+08 2e+08 3e+08 4e+08 5e+0 



2e+08 2.5e+08 
"E_32, all odd twists" + 



1e+08 2e+08 3e+08 4e+08 5e+08 6e+08 7e+08 8e+08 9e+08 
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flatness we are observing reflects the main dependence on T 3 / 4 . 
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0.5 1 1.5 2 2.5 5 10 15 20 25 30 

FIGURE 3. The first picture depicts the value distribution of 
L Ell (l/2,Xd), for prime \d\, -788299808 < d < 0, even functional 
equation, compared to Po(N,x), with N = 20. For contrast, we 
depict, in the second picture, the value distribution of L(l/2,Xd) 
(Dirichlet i-functions) for all fundamental 800000 < |d| < 1000000. 
Here, the Katz-Sarnak philosophy predicts a Unitary Symplectic 
family, and so we compare the data against Pusp(N,x), N = 5. 

In these pictures, we have renormalizcd the L values so as to 
have the same means as Po(20, x), and -FVs p (5, x) respectively, 
and have not incorporated the dk values into the pictures. 
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